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We study smooth, caustic-free, chaotic semiclassical dy- 
namics on two-dimensional phase space and find that the dy- 
namics can be approached by an iterative procedure which 
constructs an approximation to the exact long-time semiclas- 
sical propagator. Semiclassical propagation all the way to the 
Heisenberg time, where individual eigenstates are resolved, 
can be computed in polynomial time, obviating the need to 
sum over an exponentially large number of classical paths. 
At long times, the dynamics becomes quantum-like, given by 
a matrix of the same dimension as the quantum propaga- 
tor. This matrix, however, differs both from the quantum 
and the one-step semiclassical propagators, allowing for the 
study of the breakdown of the semiclassical approximation. 
The results shed light on the accuracy of the Gutzwiller trace 
formula in two dimensions, and on the source of long-time 
periodic orbit correlations. 



I. INTRODUCTION 

Semiclassical methods have a long history dating back 
to the very beginnings of quantum mechanics, and have 
provided insight into many properties of quantum me- 
chanical systems. These methods provide a bridge, ex- 
pressing quantum behavior in terms of classical paths 
and their corresponding actions. In intcgrable systems, 
the connection between quantum and classical behavior 
is well understood through EBK quantization techniques, 
which lead to an intuitive understanding of the quantum 
properties of these systems. For a nonintegrable system, 
it is not nearly as clear how much of the quantum be- 
havior (e.g. spectrum, eigenstates, long-time dynamics, 
transport) can be understood via semiclassical methods. 
One would like to be able to separate out those features 
of the quantum behavior which can thus be explained 
in terms of interference between classical paths from the 
"hard quantum effects" , such as diffraction and tunnel- 
ing. 

Although the Van-Vleck formula describing short-time 
quantum behavior in terms of classical paths has been 
around for many years, it is only starting with the work 
of Gutzwiller |jj that it has been possible to try to un- 
derstand long-time quantum properties in terms of clas- 
sical behavior. Much insight is given by the Gutzwiller 
trace formula in the energy domain, which relates the 



quantum spectrum to a sum over classical periodic or- 
bits. However, the formula is a formal expression, that 
needs to be resummed to get convergence in the limit 
where more and more orbits are included g]. Further- 
more, in practice summing over long orbits (ones with 
period comparable to Ti times the density of states of the 
system, so that individual eigenstates can be resolved) is 
not feasible, because the number of classical orbits grows 
exponentially with time in a chaotic system.^] Cycle ex- 
pansion methods, which use the symbolic dynamics of 
the underlying classical system to express long periodic 
orbits in terms of shorter ones have been very important 
in this regard ||. 

However, this still leaves open the question of the prop- 
erties of the semiclassical dynamics in the time domain. 
Much important work here has been done in systems like 
the stadium billiard and the baker's map @. Numer- 
ical evidence was produced, and theoretical arguments 
given, showing clearly that the semiclassical approxi- 
mation works well past the mixing time of the system, 
where multiple stationary paths contribute to the quan- 
tum propagator, and where the purely classical approx- 
imation (without interference effects) breaks down com- 
pletely. This was initially somewhat surprising, because 
after the mixing time the classical dynamics begins gen- 
crating structures in phase space on scales smaller than 
Planck's constant, while the quantum dynamics washes 
out any information on these scales. Thus the semiclas- 
sical approximation tries to follow the quantum prop- 
agator by summing over an exponentially large num- 
ber of paths. This exponential proliferation of classi- 
cal paths prevents one from performing meaningful long- 
time quantum-semiclassical comparisons for small values 
of h (large Heisenberg time). The fact that the quan- 
tum propagator effectively smears out all of the sub-?i 
structure in phase pace, and thus contains an amount of 
information which scales only as a power law in fi sug- 
gests that some such reduction may also be possible in 
the semiclassical calculation. This is important because 
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1 Special cases are known for which the long-time semiclas- 
sical dynamics can be computed exactly, e.g. the cat map 
and geodesic flows on surfaces of constant negative curva- 
ture. However, the semiclassical evolution for these systems 
is in fact the same as the quantum evolution because of their 
linearity (or homogeneity in the case of constant negative cur- 
vature), making them less interesting as a testing ground for 
the general applicability of semiclassical methods. 
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a priori it is not at all obvious (1) that the long-time 
semiclassical dynamics converges in the sense of having 
well defined stationary states at the Heisenberg time, or 
for that matter (2) that in those cases where the semi- 
classical dynamics does converge, what it converges to 
actually approximates the quantum dynamics, station- 
ary states, etc. The breakdown of semiclassical validity, 
in particular in the presence of diffraction, discontinu- 
ities, and caustics, is also a somewhat controversial issue 
which has been difficult to address in practice because of 
the computational obstacles. 

A key reason for the finiteness of long-time informa- 
tion in the quantum case is the multiplicativity of the 
quantum propagator. In other words, the long-time dy- 
namics can be obtained by simply iterating the short- 
time propagation (matrix multiplication in the case of a 
finite-dimensional Hilbert space). The semiclassical dy- 
namics does not share this property because concatena- 
tions of classical paths in general produce paths that are 
not classical. Stationary phase integrals must be per- 
formed to obtain the long-time semiclassical propagator 
from shorter steps, and the number of stationary phase 
points grows exponentially with time. An interesting 
way to make the semiclassical dynamics multiplicative 
by extending the space on which the multiplicative oper- 
ator acts has been described by Cvitanovic and Vattay 
M. It has also been seen that in the special case of the 
baker's map, the long-time dynamics can be computed 
with good accuracy in polynomial time using the Heisen- 
berg uncertainty principle and the exponential decay of 
time correlations in chaotic systems || 7| . These methods 
(of which three are known to the author) effectively col- 
lect together all contributions from classical paths which 
come together on scales much smaller than h. Semiclas- 
sical amplitude thus collected can then be propagated 
again, thus making the resulting dynamics multiplica- 
tive. All the methods, however, made use of the simple 
symbolic dynamics of the baker's map and its very spe- 
cial structure in position and momentum space. Some 
of the ideas of consolidation on sub-?i scales, however, 
may yet turn out to be fruitful in analyzing more generic 
systems, with no such special structure 0. 

In this paper, we take a somewhat different approach 
to the problem, inspired by the fact that in the baker's 
map it was found pj that a good approximation to the 
long-time semiclassical dynamics could be obtained even 
if one periodically projects the higher-dimensional semi- 
classical vector onto the iV-dimensional quantum space. 
In other words, the true semiclassical dynamics is eval- 
uated exactly for Tq steps. Matrix elements are then 



taken between quantum states, producing a matrix of the 
same dimension as the quantum propagator. This ma- 
trix is then iterated to produce an approximation to the 
long-time semiclassical behavior. Now for Tq — 1 this is 
very reminiscent of Bogomolny's approach to quantiza- 
tion ||. It should be noted, however, that for Tq = 1, 
this procedure does not produce anything at all resem- 
bling the long-time semiclassical dynamics. In fact for 
small Tq one gets closer and closer to the exact quantum 
dynamics, as can be seen by taking the limit Tq — > and 
recovering the Feynman form of the quantum propaga- 
tor. This is good if one just wants to know the quantum 
answer, but not if one is interested in understanding how 
quantum dynamics is affected by classical properties, in 
the sense, for example, of the Gutzwiller trace formula. 
What is interesting here is that for Tq > 1, but still 
much smaller than the Heisenberg time, one can obtain a 
good approximation to the actual long-time semiclassical 
dynamics, and thus study the validity (and breakdown) 
of the semiclassical propagation. Detailed comparisons 
become easily possible between quantum and semiclassi- 
cal spectra, eigenstates, long-time transport, and other 
properties. 

Here we apply these ideas to a somewhat different class 
of systems, ones that lack the somewhat unpleasant "cut 
and paste" discontinuities of the baker's map, that are 
not piecewise linear, and which are also more generic 
in the sense of constituting a large family and allow- 
ing for continuous perturbations. These systems, namely 
smooth automorphisms of the torus (e.g. kicked maps, 
perturbed cat maps) lack a special symbolic dynamics 
and have no preferred basis, making them a good test 
case for the ideas outlined abovej^] It should be pointed 
out that the procedure presented in what follows does 
not make use of detailed knowledge about the system at 
hand. In this way, it may be more robust than, for exam- 
ple, the cycle expansion methods in periodic orbit theory. 
Because no explicit use is made of periodic orbits, it need 
not be the case here that long periodic orbits can be ex- 
pressed in terms of short periodic orbits, for example. All 
that is needed is that long paths can be constructed out 
of short paths, something that is true in all but the most 
pathological cases. 

The rest of this paper is organized as follows: we first 
briefly review the classical, quantum, and semiclassical 
dynamics of the systems under consideration. The iter- 
ation of semiclassical propagators is discussed, and error 
estimates are given as a function of propagation time, 



2 Tq is called the quantization time, and can be thought of 
the time domain analogue of the T*parameter in the spectral 
theory of Bogomolny and Keating |8|. One should note, how- 
ever, that in periodic orbit theory, the critical time T* should 



scale as the mixing time, logarithmically with Ti, whereas our 
scale Tq is of order one in units of the shortest periodic orbit. 
See the discussion in Section^. 

3 An important constraint is that we shall be looking exclu- 
sively here at caustic free, purely chaotic systems. Caustics, 
discontinuities, and mixed phase space are all important sub- 
jects of inquiry and will be looked at in a future paper. 
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quantization time, and h, with particular emphasis on 
Heisenberg-time propagation. Numerical tests follow, 
showing the convergence of the iterative approximation 
as well as measuring the deviation from exact quantum 
and one-step iterative results. Then the "effective one- 
step semiclassical propagator" is discussed, which with 
exponential accuracy describes long-time semiclassical 
propagation, and which differs from the quantum and 
one-step semiclassical matrices. Again, numerical find- 
ings are presented. Finally, the conclusion treats some 
general questions and addresses possible applications and 
extensions of the results obtained. 

In a companion paper to this article |l0| ], the meth- 
ods are extended and applied to the study of semiclassi- 
cal dynamical localization in classically diffusive systems, 
showing that in fact interference between classical paths 
is sufficient to understand the end of quantum diffusion 
in these systems. 



II. THEORETICAL ANALYSIS 

A. Hard chaotic dynamics on a torus 

We will consider smooth, chaotic classical dynamics on 
a compact phase space, satisfying the Anosov property,^] 
and free from caustics in the coordinate system of inter- 
est. To be specific, let us take the following family of 
discrete area-preserving maps on a torus: 



p — > p = p + mq — V (q) mod 1 
q — > q = q + np + T (p) mod 1 . 



(1) 



The above dynamics can be obtained from the strobo- 
scope discretization of a kicked system with a kick 
potential — \mq 2 + V(q) applied once every time step and 
a free evolution governed by the kinetic term ^np 2 +T(p). 
Here m, n are arbitrary integers, while V, T axe smooth 
periodic functions. The system can also be thought of as 
a continuous perturbation of the linear system (cat map) 

M 



p — ► p = p + mq mod I 

q — ► q = np + (mn + l)q mod 1 . 



(2) 



The Jacobian of the transformation in Eq. |l] is given 



by 



J 



1 m-V"(q) 
n + T"{p) l + ( n + T"(p))(m-V"(q)) 



(3) 



We notice that for given integers m, n we can choose the 
functions V and T such that the quantity m — V" (q) is 
strictly greater than for all q and similarly for n+T"(p). 
(In words, we ensure that the system everywhere looks 
locally like an inverted harmonic oscillator.) This im- 
plies strict positivity of all four entries in the Jacobian 
matrix, a property that is of course preserved under itera- 
tion of the dynamics. Furthermore, positivity of both off- 
diagonal entries implies hyperbolicity, because it ensures 
det J > 2. All such systems therefore provide examples of 
hard chaos, being free of integrable regions. To see that 
they are free of caustics in either position or momentum 
space, it is sufficient to note that the two off-diagonal en- 
tries of the Jacobian, dp(t)/dq(0) and dq(t) / dp(0) always 
remain non-zero. 

The quantization of kicked systems is straightforward 
and well-covered in the literature [[flj . We take % so that 
N = 1/27t7j has an integer value. Then an TV-dimensional 
position basis for the Hilbert space is given by \qi), where 
qi = (i + 6q)/N, i = 0...N — 1. Similarly, the mo- 
mentum space basis is given by \pj), with allowed values 
Pj = (j + Ei)/N, j — 0. . . N — 1. e ,i form a family 
of possible quantization conditions (they correspond to 
phases associated with circling the torus in the p and q 
directions, respectively). The two bases are related by a 
discrete fourier transform. The dynamics is now defined 
by the unitary N x N matrix 



U = cxp 



exp 



i ( \np 2 + T(p) ) /Ti 



- i [ hnq 2 -V(q) ) /h 



(4) 



where each factor is evaluated in the appropriate basis, 
and an implicit forward and backward fourier transform 
has been performed.^ The key point is that U here is just 
a matrix, and the long-time quantum evolution of the 
system is given by matrix multiplication. Alternatively, 
the matrix U can be diagonalized to find the stationary 
properties (eigenstates and eigenvalues) of the system. 
All this can be performed in polynomial time, and the 
total amount of information contained in the quantum 
system is of order N 2 . In particular, all of the informa- 
tion is present in the one-step quantum propagator. 



B. Semiclassical dynamics 

We now consider the long time semiclassical dynam- 
ics as given by the Gutzwiller-Van-Vleck propagator jl] 
evaluated between quantum mechanically allowed states. 
In position space, the propagator has the form 



4 This means that at each phase space point the tangent 
space can be decomposed into an exponentially expanding 
linear subspace and an exponentially contracting linear sub- 
space. For two-dimensional phase spaces, there is one expand- 
ing and one contracting direction at each point. 



5 As in the linear (cat map) case, continuity of the potential 
and kinetic terms (mod Ti) requires that iV be chosen even, 
unless both m and n are even. 
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G sc (q,q',t) 



2nih 



d/2 



E 



1/2 



x exp 



iSj(q,q' ,t) iirvj 



(5) 



where Sj is the action for classical path j, the determi- 
nant corresponds to the classical probability density for 
going between q and q' via this path, and the phase is 
given by the action, corrected by the count of conjugate 
points Vj. For any given time t, this produces an JV x JV 
matrix A t . The semiclassical evolution must of course be 
evaluated in some coordinate system, such as position or 
momentum, but once the matrix has been constructed 
it can be rotated into any quantum basis that one finds 
convenient. Being a matrix connecting quantum in and 
out states, A t looks like a quantum object, but it is not 
unitary (A t A- t 7^ I) and does not satisfy multiplicativ- 
ity {At+t> 7^ A t At'). Thus it is not a priori obvious, for 
example, that the semiclassical evolution at long times 
converges to a well-defined set of eigenstates and eigen- 
values, i.e. A t ip n ~ s^ipn, or whether such eigenvalues 
and eigenstates e n , ip n have any connection with those of 
the quantum dynamics. 

Let us study the deviation from multiplicativity of the 
semiclassical dynamics in a smooth chaotic system. We 
note first that for semiclassical propagators in the time 
domain, 



A 2t (q",q)= / dq'A t (q",q')A t (q',q), 



(6) 



sp 



where J indicates that the equality holds only if the in- 
termediate integration is performed by stationary phase 
fl3f . If the integration is performed exactly instead of in 
the stationary phase approximation, we obtain a relative 
error of order Ti in the answer.^ This error comes from 
higher-order terms in the stationary phase expansion. (In 
the presence of caustics, the coefficient of the 0(K) term 
can blow up in certain regions of space, eventually dom- 
inating the semiclassical evolution. Possibilities for han- 
dling this problem include uniformization or choosing a 
different basis (e.g. a gaussian basis) for performing the 
semiclassical calculation. To avoid these serious difficul- 
ties we will deal throughout this paper with systems that 
are caustic-free in the chosen basis.) 

For systems on a compact classical phase space, po- 
sition and momentum values are of course labeled by 
discrete integers, so strictly speaking the notion of sta- 
tionary phase integration is not well defined. To make 



6 Of course, for large t, many stationary paths are summed 
over on the right hand side of Eq. ^. The relative error in 
ignoring subleading terms in the stationary phase expansion 
for each such path is 0(h). So the fractional error in the full 
answer is also 0(h), provided that the errors add no more 
coherently than the leading terms themselves. 



sense of the semiclassical dynamics, one must rewrite the 
sums over topological classes (winding numbers) using 
the Poisson summation formula, and evaluate the result- 
ing integrals by stationary phase |l2]. This produces a 
semiclassical propagator defined only on a discrete po- 
sition (or momentum) grid. However, in the limit of 
small Ti (large N), the spacing between quantum basis 
states vanishes, and the discreteness of the quantum ba- 
sis ceases to be physically significant. This is true as long 
as all quantum structures in phase space are "generic" , 
scaling as y/h in both the q and p directions, while the 
spacings Aq and Ap scale as Ti. The error analogous to 
Eq. I that we are interested in for a compact phase space 
is the difference between performing a sum over inter- 
mediate channels q' and performing the corresponding 
integral by stationary phase. We can then write 



JV-l 



A 2 S, *0 = E Mi,j)Mi, k) + 0(N-^ 2 ) . 



(7) 



3=0 



Note that the error scales as N~ 3 ^ 2 ~ fi 3 ^ 2 because nor- 
malization (probability conservation) requires that the 
actual matrix elements of A t and Ait scale as N~ x / 2 ~ 
h 1 / 2 . The relative error is scaling as TV" 1 ~ h. We now 
define a natural norm for measuring the difference be- 
tween two matrices (under which a unitary matrix has 
norm one), 



\A-B\ 



L tx{A -B)\A- 



B) 



(8) 



We then have 



\Ao 



(At) 



2 1 1 2 



o(h 2 ), 



(9) 



whereas the matrices A t individually have norm of order 
unity. Of course, the coefficient in front of the 0(h 2 ) 
depends on the amount of nonlincarity in the underly- 
ing classical dynamics. For example, the cat map given 
by Eq. |^ is exactly linear, and produces a semiclassi- 
cal dynamics A t which is exactly multiplicative. This, 
however, is not very interesting because in that case A t 
is also equal to the exact quantum dynamics 17*. For 
generic perturbing potentials V{q) or T(p), the semiclas- 
sical answer differs form the quantum, and in that case 
the number multiplying h 2 will indeed be of order unity. 

The result in Eq. ^ is already quite promising. It tells 
us that for small %, it is a very good approximation to 
compute the semiclassical dynamics exactly for 10 steps 
and then square the matrix instead of trying to do the 
exact calculation for 20 steps. The former is a much 
easier problem to solve because the number of classical 
paths needed to compute A t scales exponentially with t. 
Inspired by this, we ask to what extent we may approx- 
imate the exact time-i semiclassical propagator for large 
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t by dividing t into more and more shorter time inter- 
vals. Replacing A t by [A t i M ) M involves making an error 
M—l times, each time approximating a stationary phase 
integration by exact multiplication. Assuming these er- 
rors add incoherently, and taking M to be large, we then 
have 

\\A t -(A t/M ) M \\ 2 =0(Mh 2 ). (10) 

One can also obtain this result by iterating the procedure 
indicated by Eq. ^| for M that is a power of two. (In other 
words, we express A t in terms of A t / 2 plus an error term, 
then A t /2 in terms of A t /^ plus another error term, etc.) 
However we need to be careful about the assumption of 
incoherent accumulation of errors. Because the Hilbert 
space is finite-dimensional, eventually we must consider 
interference between different error terms. To see this in 
our formalism, let B = A t/M and C = {A t ) 1/M . (Al- 
though the M-th root of a matrix is in general an am- 
biguous quantity, here there is no ambiguity in what we 
mean by the matrix C : we simply choose that root which 
is closest to B.) Now the quantity that we are interested 
in is || C M - || 2 . Let e = C - B. To lowest order in e, 

||C M - B M \\ 2 = \\eB M ~ 1 + BeB M - 2 + ... + B M - 1 e 
+ 0(e 2 )\\ 2 . (11) 

Now we work in the basis in which B is diagonal, and 
write e as the sum of its diagonal and off-diagonal parts 
in that basis. The diagonal part of e commutes through 
B, giving a coherent contribution from the M terms in 
the sum. The off-diagonal part leads to an incoherent 
contribution because the eigenphases of B are generic. 
So we obtain for the basis-independent norm 

\\C M - B M \\ 2 = 0(A/||e off _ diag || 2 ) + 0(M 2 ||e diag || 2 ) 

+ 0(Af 2 ||e 2 || 2 ) + ... . (12) 

Now as long as the errors in the stationary phase approx- 
imation are not preferentially diagonal (i.e. not in gen- 
eral tending to multiply the exact answer), the weight 
of the matrix e that is on the diagonal is a fraction 
1/N ~ h of the total weight of e. Noting from Eq. ^ that 
||e|| 2 = 0{h 2 ), we see that ||e diag || 2 = 0(h 3 ). Further- 
more, for the values of M that we are going to consider 
(up to the Heisenberg time Tji ~ Ti _1 ), terms higher 
order in e (e.g. 0(M 2 Ti )) can be ignored. So we obtain 

\\A t - (A t/M ) M \\ 2 = 0(Mh 2 ) + 0(M 2 h 3 ) , (13) 
valid for M < 0{h~ l ). 

C. Heisenberg time dynamics 

In particular let us consider what happens at the 
Heisenberg time t = Th ~ 0{h~ l ) (computing the semi- 
classical dynamics to times longer than this will not 



produce interesting new information because individual 
eigenstates and eigenvalues will already have been re- 
solved). Let Tq be the "quantization time", the time for 
which we will compute the semiclassical propagator At„ 
exactly. First, letting M = t/TQ, we can rewrite Eq. [l3| 
as 

Then by the Heisenberg time t = O^" 1 ), we accumulte 
an error 

Pt h -(A Tq )^|| 2 = o(A^ +0 ^ . (15 ) 

The above formula is expected to hold for all Tq > 1, 
where classical paths exist connecting any two coordi- 
nate points. If we apply the result to the case Tq = 1, 
which characterizes the one-step Bogomolny propagator, 
we find 

\\A TH -{A l f»f = 0{n). (16) 

So one result of the above calculation is that for smooth 
one-dimensional kicked systems (or for two-dimensional 
Hamiltonian systems, where the scaling of the density of 
states works identically), in the absence of caustics the 
exact semiclassical propagator deviates from the iteration 
of the one-step dynamics only by an error term of order 
Vh, by the Heisenberg time. How does this propagation 
relate to the exact quantum dynamics? For a single time 
step we know that the relative error between semiclassical 
and quantum amplitudes scales as h, so 

\\A 1 -U\\ 2 = 0(h 2 ). (17) 

Following a line of reasoning completely analogous to the 
one that took us from Eq. |] to Eq. [l6] (i.e. noticing that 
off-diagonal terms in the error add incoherently, etc.), we 
find 

||(Ai) T *-[/ T *|| 2 = 0(7L). (18) 

Now combining this with Eq. [l6] we obtain a relation 
between the exact semiclassical and exact quantum an- 
swers: 

\\A TH -U T »\\ 2 = 0{h). (19) 

Thus, in the absence of discontinuities and caustics the 
semiclassical approximation is expected to do quite a 
good job in approximating the quantum dynamics all 
the way out to the Heisenberg time (in two dimensions, 
that is - in three dimensions the same analysis leads 
to the conclusion that the semiclassical approximation 
is marginal at the Heisenberg time). 
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But there still is a difference between the long-time 
semiclassical and quantum behavior, implying also a dif- 
ference between the corresponding eigenvalues and eigen- 
states. These corrections arise from "hard quantum" ef- 
fects, those beyond the stationary phase approximation. 
To be able to separate these from those effects which are 
purely semiclassical one needs to be able to compute At 
explicitly for t ~ Th, to a better approximation than 
that given by the exact quantum mechanics, Eq. [l9| It- 
eration of the one-step semiclassical propagator A\ will 
not do, since it differs from the true semiclassical propa- 
gation as much as from the quantum. So we proceed to 
consider a quantization time Tq > 1. From Eq. |l5| we 
see that iteration of At q comes much closer to the exact 
semiclassics than the difference between the latter and 
quantum mechanics. 



\t/T Q 



A* 



« \\A t - W 



(20) 



for t Th- This implies 

\\(A Tq Y^ - A t f « \\{A TQ fl T « C/*|| 2 , (21) 

and thus iteration of At q gives us an answer much closer 
to At than to 17*, allowing us to see how the long-time 
semiclassical dynamics differs from the quantum dynam- 
ics, without the need for doing an exponentially large 
amount of work. The approximation is a controlled one, 
and one can keep increasing Tq until the desired level of 
convergence to the exact semiclassics is reached. 

These ideas, extended properly to quantization of non- 
compact phase spaces, can be used for example to see 
dynamic semiclassical localization without having to sum 
explicitly over a number of classical paths that is expo- 
nential in the localization time. The results are presented 
in a companion paper [ jTo| , and address a long-standing 
question in the literature over whether dynamic localiza- 
tion is a semiclassical or hard quantum phenomenon. 

III. NUMERICAL TESTS 

We proceed now to justify numerically the power- 
counting arguments presented in the previous section. 
The system we will use is the kicked system of Eq. [|, with 
m = n = 1, and V(q) — — (K/(2tt) 2 ) sin 2irq. This can 
be thought of as a standard map (kicked rotor) with an 
extra inverted harmonic oscillator potential 1 

a sinusoidal perturbation of the 1 



1 2 



cat map. T(p), 



the perturbation of the quadratic kinetic term, is set to 
zero for simplicity. Also for simplicity periodic boundary 
conditions, with no phases, are imposed, i.e. eo = e± = 0. 
The caustic- free condition requires kick strength \K\ < 1. 
The system is then guaranteed to satisfy the Anosov 
property, as explained in the discussion following Eq. |^. 
We choose N = 256 to be the dimension of the Hilbert 
space. This is well in the semiclassical regime, and large 



enough so that a direct summation of classical paths to 
the Heisenberg time (there are 3 256 of them) is clearly 
not practical. 

We then fix a value of K, at K = 0.5, and compute the 
matrix elements of the semiclassical propagator between 
quantum states in the momentum basis (using p is nat- 
ural because we are thinking of this as a kicked system; 
also working in momentum space leads naturally into 
the study of dynamical localization). Exact semiclas- 
sical matrices Aj- Q are computed for quantization times 
Tq = 1 . . . 8, using the Gutzwiller-Van-Vleck expression, 
Eq. ^. The convergence of the iterative approximation 
can then be investigated. We first compute the quan- 
tity ||(A r? )*/ T « - (A 8 )*/ 8 || 2 for T Q = 1,3,5,7, and for a 
range of times t that extends beyond the Heisenberg time 
Th = 256. As explained above, the exact semiclassical 
propagator for times of order Tjj cannot be computed 
exactly^ , so we use the difference between approxima- 
tions at different values of Tq as measure of convergence 
as Tq — > oo. We assume errors in successive approxima- 
tions are uncorrelated, i.e. 
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Tq 
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t/T'c 



Q II 2 
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Tq 
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1-4 



-M 



(22) 



For comparison purposes, the difference with the quan- 
tum dynamics \\Ut — (^4s)'^ 8 || 2 is also computed. All these 
are plotted as a function of time t in Figure [l]. From 
top to bottom, the quantities plotted are the squared 
differences between the eight-step iterated semiclassics 
A 1 ^ 8 and (i) the one-step iterated semiclassics A\, (ii) 
the quantum mechanics 17* , (iii) the three-step iterated 

semiclassics A^ 3 , (iv) the five-step iterated semiclassics, 
and finally (v) the seven-step iterated semiclassics. Each 
set of points is also fitted to a function of the form at+bt 2 , 
as suggested by Eq. [l4|. 

We notice first of all that the three-step, five-step, and 
seven-step approximation come progressively closer to 
the eight-step approximation which is our basis of com- 
parison. The differences between all these are signifi- 
cantly smaller than that between any of these and the 
quantum dynamics. Finally, the one-step iterated ap- 
proximation does a (relatively) poor job of reproducing 
the long-time semiclassical behavior (however the differ- 
ence between it and the other calculations is still small, 
due to the smallness of ft) . Consistent with the prediction 
of Eqs. [l4], |2^, all the curves are well fitted by the sum of 
a linear and quadratic functions of time. Furthermore, 
at least for Tq > 3, the linear term is seen to dominate 



7 At the end of the next section, we will see an example of a 
numerical test for moderate time t, where exact semiclassical 
calculation is in fact possible, though very time-consuming. 
There, approximation techniques are explicitly shown to pro- 
duce a very good answer with much less effort. 
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for times up to the Heisenberg time (Tjj — N = 256 in 
this calculation). Thus, to understand the convergence 
of the iterative method for Tq ^> 1 , it is sufficient to look 
at the behavior of the linear term in Fig. [3]. 

This is in fact done in Fig. B. For each Tq , we find the 
numerical value of cl(Tq) that fits 



I Tq 



Al 



a(T Q )t + b(T Q )t 2 



(23) 



In Fig. ||, the quantity cl{Tq) is plotted vs. Tq, using 
plusses for N = 256 and crosses for N = 128. However, 
for Tq = 1, the difference between the quantum dynamics 
and the eight-step iterated semiclassics is plotted, i.e the 



linear coefficient of ||[7* — Ag^ 8 || 2 , not \\A 



Al /S \\ 2 . The 

latter value is off the scale in the Figure, being 30 • 1CT 7 
for N = 256 and 115 • 1(T 7 for N = 128 (notice that this 
error, though relatively large compared to those obtained 
for bigger Tq, still has the right h 2 dependence). From 
Eqs. [LJ, 03, we predict 



1 



1 



a{T Q )=ch\- + —) 

° J-Q 



(24) 



for Tq 3> 1, where c is an undetermined constant of order 
unity. This predicted behavior fits the observed values of 
q,(Tq) quite nicely, with c = 0.29. The resulting curves, 
following Eq. ^4| are plotted in Fig. ||. We see that the er- 
ror in the iterative approximation indeed scales as Ti 2 per 
time step, leading to an error of order Ti by the Heisenberg 
time. Furthermore, the error goes to zero as Tq 3> 1, in 
the predicted manner, and is already much smaller than 
the difference between the quantum mechanics and the 
semiclassics by Tq = 3. 



IV. EFFECTIVE ONE-STEP LONG-TIME 
PROPAGATOR 

We have seen in the preceding sections how to obtain 
convergent approximations to the semiclassical propaga- 
tor for long times, including times beyond the Heisen- 
berg time. This dynamics can now be fourier trans- 
formed to obtain local densities of states for various 
initial wavepackets, and from these the eigenstates and 
eigenvalues could be extracted. In particular, the semi- 
classical spectrum can be obtained by tracing over the 
fourier transformed semiclassical evolution. The quanti- 
zation time Tq can be increased until the desired level 
of convergence is reached, and the eigenstates and eigen- 
values obtained in this way can finally be compared with 
those of the quantum matrix. 

This, however, is a rather tedious process, requiring the 
approximate semiclassical propagator to be evaluated at 
many values of t, from short times to times well beyond 
Tjj- Moreover, if we settle on a fixed quantization time 
Tq, and computed all long-time dynamics using it, we 
will finally obtain nothing more than the eigenstates and 
eigenvalues of the matrix At q . So clearly the sensible 



thing to do is to diagonalize At q directly for a range 
of values of Tq and check for convergence as Tq ^> 1. 
Unfortunately, a technical difficulty arises here. Eigen- 
states of the true dynamics (quantum or semiclassical) 
with eigenphases separated approximately by an integer 
multiple of 2tt/Tq may get mixed with each other in the 
diagonalization of At q , preventing one from extracting 
the true eigenstates. This behavior appears to be generic 
(it does not happen in A\ due to level repulsion in the 
presence of chaos). One can try to get around this dif- 
ficulty by comparing the eigenstates of At q for a range 
of Tq, and selecting those that come closest to agreeing 
for several values of Tq. Such a procedure does in fact 
appear to converge to a reasonable set of N eigenstates, 
and it might be expected to produce results comparable 
to those that would be obtained by fourier transforming 
long-time dynamics produced by using several matrices 

However, a simpler solution now presents itself. Given 
that we believe At q have similar eigenstates for all values 
Tq ^> 1, these must also be the eigenstates of 



D Tq = A Tq A Tq _ 



(25) 



for example. No mixing of eigenstates should arise in 
Dt q (because of level repulsion), and moreover this 
method has the advantage that the semiclassical eigen- 
values can be read off directly, without having to decide 
which root of an eigenphase must be taken in each case. 
Tq is then increased, and convergence to the "true" semi- 
classical eigenstates and eigenvalues is obtained (conver- 
gence is expected based on what we know about the con- 
vergence of the dynamics from the two previous sections). 

In fact, it is sufficient to look at convergence of the 
Dt q matrix. In Fig. the results of such an analysis are 
presented, for t he same system as that studied numeri- 
cally in Section III . Here a Heisenberg time (matrix size) 



N = 64 is used. We plot the squared norm \\Dt q — Dg\\ 
for 2 < Tq < 6 (for Tq > 6, numerical errors begin 
to play a role in the error analysis, however, the down- 
ward trend continues to 5.4 • 10~ 12 for the difference be- 
tween D$ and Dg). For comparison, we also measure 
Pi - D 9 \\ 2 = 4.6 ■ 10- 5 and \\U - D 9 \\ 2 = 6.2 ■ lO" 6 . 

Remarkably, we find exponentially fast convergence 
with the quantization time, in marked contrast with the 
power law convergence obtained previously for the dy- 
namics. The data agrees well with the exponential form 
h 2 exp(8.1 - 3.2T Q ), which is also plotted in Fig. [| This 
exponential behavior is consistent with analogous results 
seen in cycle expansion methods. How can it be recon- 
ciled with the power-law behavior seen in Eq. [13? Let 



8 For example, we could approximate Aioo ~ AfAj, Aioi ~ 
-AgAg, etc. This procedure is in fact useful for evaluating 
long-time propagators at arbitrary times t, including those 
that are not divisible by a suitable value of Tq. 
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£>*= lim D Tq . (26) 

Tq — »oo 

Now clearly A t = D*(l + e*), where et is a correction 
falling off exponentially with t, would satisfy the find- 
ing (upon substitution into Eq. |2^) that Dt q converges 
exponentially quickly. However, it would not be consis- 
tent with Eq. [L|, which requires power-law behavior for 
the dynamics, and which we have seen verified numeri- 
cally. The two results can be reconciled if we notice that 
exponential convergence of Dt q to D* only requires 

A t = Dl{l + e* + e t ). (27) 

Here ||et|| 2 = 0(h 2 exp (— at)), while e* is a 
t— independent matrix with norm He*!) 2 = 0(h 2 ). a 
is a constant associated with the mixing time scale of 
the underlying classical dynamics. Now when taking the 
product D t = A t At\, the e* contribution cancels, pro- 
ducing D+ plus a time-dependent correction that goes as 
0(K 2 exp (— at)). However, if we compare A t with A^j M 
using Eq. where t and t/M are assumed for simplic- 
ity to be large, we notice that the latter quantity has an 
extra M — 1 terms of order compared with the exact 
semiclassics At- This is consistent with what we found 
in Eq. [If], and leads to power-law errors as in Eq. [l^. 
For example, from Eq. |2^ we see Ait — A 2 — 0(e*), in 
agreement with Eq. [)| 

Our interpretation of Eq. ^7] is that e t is an error that 
results from concatenating short classical paths to pro- 
duce long ones. It is thus analogous to the errors obtained 
in approximating long periodic orbits by short ones in cy- 
cle expansion methods, e*, on the other hand, is an error 
in some way associated with projecting the full scmiclas- 
sical dynamics for any time t onto quantum initial and 
final states. This error is therefore independent of the 
time after which such projection is performed. It should 
be possible to make a connection between this effect and 
what happens in the baker's map, for example, when a 
slightly higher-dimensional effective semiclassical space 
(in which the semiclassical dynamics is almost exactly 
iterative) must be projected onto N— vectors to produce 
quantities that can be compared with quantum matrix 
elements ||. One would also like to understand better 
the parametrization ambiguities in the definition of D+ 
and e*. For example, we could have written 

A t = (1 + e* + l t )Dl (28) 

as an alternative to Eq. |27|. This would be natu- 
ral is we had looked at the quantity A^^Atq (= 

(A-Tq^i-TqV) instead of Dt q as defined above. 

In any case, we see that the matrix D+ is key to un- 
derstanding the long-time semiclassical dynamics. This 
matrix is an effective one-step propagator that can be 
used to obtain the semiclassical propagator at time t + 1 
with exponential accuracy, given the propagator for time 
t. It thus gives to us in a trivial way the stationary prop- 
erties of the long-time semiclassical evolution. 



To show how the factorization given in Eq. |27j can be 
used in practice, we take the kicked map studied earlier 
in this section and compute an approximation to Ag by 
evaluating 

AgxiDKl + e*) 
= £>^(l + e*) 
*DlA 5 

w (A 5 A^) 4 A 5 (29) 

All errors e t in Eq. ^9| are exponentially small in the rel- 
evant times t: an error of size eg is made in the first line, 
and errors of size £5, £4 in the third and fourth lines, the 
latter of course dominating the error in the final answer. 
Notice that we need to evaluate the semiclassical dynam- 
ics A t exactly at two values of t to obtain the two matrices 
and e*. After computing the semiclassical propaga- 
tor Ag exactly, we find (factoring out the h dependance), 
that the error |j {A 5 A^ 1 ) 4 A 5 - A>|| 2 is 0.00003?i 2 . In the 
same units, we find the error HA5A4 — Ag\\ 2 to be 0.5?i 2 , 
still not bad, but lacking the exponential error suppres- 
sion. For comparison, the difference between the exact 
semiclassical 9— step propagator Ag and the quantum me- 
chanics Ug is 3. IS 2 , while the difference from the 9— times 
iterated one-step propagator {\\A\ — Ag\\ 2 ) is 89.5?i 2 . (It 
is interesting to note just how poorly iteration of the one- 
step semiclassical propagator A\ docs in reproducing the 
correct longer-time semiclassical behavior.) 

Of course, this procedure can now be extended to times 
much longer than 9, where exact evaluation of the semi- 
classics is clearly not practical. Six digits of accuracy in 
long-time semiclassical propagation (beyond what is ex- 
pected merely from the smallness of K) are obtained with 
only a modest amount of work. 

V. CONCLUSION 

We have seen in the preceding section that long-time 
semiclassical behavior is given to exponentially good ac- 
curacy by iteration of an effective matrix D± that is dif- 
ferent from both the quantum evolution matrix U and 
the one-step semiclassical propagator A\. In fact, the 
one-step semiclassical propagator contains essentially no 
information about long-time semiclassical dynamics, ex- 
cept of course for the similarity based on both being re- 
lated to the quantum evolution. However, knowing the 
semiclassical dynamics for times Tq 3> 1 (measured in 
units of the shortest periodic orbit time), enables one to 
deduce stationary semiclassical behavior, including dy- 
namical information to the Hciscnberg time and beyond. 
Why is long-time semiclassical information exhausted at 
Tq 3> 1? In periodic orbit methods, one expects knowl- 
edge of long-time dynamics to be contained in orbits of 
length up to the mixing time (~ log AT), at which point 
phase space has been explored at the scale of Planck's 
constant. Longer orbits can be produced from these if we 
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allow smearing over sub-planck coordinates. The O(N) 
periodic points at the mixing time contain 0(N 2 ) pieces 
of information if we record where each lies in relation to 
some quantum basis. Similarly, the full quantum theory 
can of course be described by an N x N matrix. The 
information required to obtain the semiclassical theory 
is only slightly bigger, scaling in the same way with TV. 
We note that collecting and iterating semiclassical infor- 
mation after one time step on a scale of 1/N allows one 
to resolve periodic orbits up to the mixing time, while 
collecting after a slightly longer exact evolution is equiv- 
alent to having knowledge of periodic orbits longer than 
the mixing time, where exponential convergence is ex- 
pected. 

The results obtained here may shed light on the ac- 
curacy of the Gutzwiller trace formula in the energy 
domain. Though the relationship between the two ap- 
proaches is a nontrivial one (the trace formula focusing 
on periodic orbits, and thus corresponding to a station- 
ary phase integration of the dynamics), one may hope 
that further progress may be made in bringing together 
dynamical and periodic orbit methods. In any case, we 
have already seen that the stationary phase integration 
implicit in the trace formula is not necessary for obtaining 
sensible semiclassical eigenvalues and eigenstates, which 
compare well with those of an exact quantum calculation. 

More work is also needed in bringing together the ap- 
proximation methods discussed here with other semiclas- 
sical consolidation techniques, such as those that have 
been used quite successfully for the baker's map HQ. 

These iterative methods may in addition provide a new 
perspective on the long-standing questions of classical 
periodic orbit correlations and of how semiclassical dy- 
namics can reproduce Heisenberg-time quantum behav- 
ior @. All such issues become much less mysterious 
once we realize that long orbits are given to a good ap- 
proximation by concatenating shorter classical trajecto- 
ries. Thus, long-time semiclassical evolution can produce 
delta-function peaks in the spectrum for essentially the 
same reason as the quantum mechanical evolution, i.e. 
it is given (approximately) by iteration of an (almost) 
unitary finite matrix. 

Many of the findings may be modified in the presence 
of discontinuities, strong diffraction, and caustics, where 
classical-quantum correspondence is more poorly under- 
stood. In particular, the proliferation of caustics may 
be expected to cause a divergence between semiclassical 
and quantum dynamics well before the Heisenberg time 
is reached M. In these situations, a better understand- 
ing of iterative dynamical methods may still allow one to 
follow the semiclassical dynamics past the scale at which 
this divergence occurs, and to see explicitly when and in 
what way the semiclassical approximation breaks down 
(and what corrections to the semiclassical formulas may 
be necessary to restore long-time correspondence). In 
some cases, certain qualitative features of the quantum 
dynamics (e.g. dynamical localization |lCj[ ], level spacing 
statistics) may be well reproduced by the long-time semi- 



classical dynamics, even when detailed correspondence 
between the propagators has been lost. 
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FIG. 1. Convergence to long-time semiclassical behavior 
with increasing quantization time Tq. A long-time calculation 
with Tq = 8 is used as the reference. From top to bottom, the 
five sets of data represent the squared difference between this 
calculation and (i) the iterated one-step semiclassics, (ii) the 
quantum mechanics, (iii) the three-step iterated semiclassics, 
(iv) the five-step iterated semiclassics, and (v) the seven-step 
iterated semiclassics. The Heisenberg time is Th = 256. Each 
curve is fit to the sum of a linear and quadratic function of 
time, in accordance with Eq. |lj. 

FIG. 2. Convergence of the linear coefficients in the previ- 
ous Figure (error per time step), with increasing quantization 
time Tq. From left to right, the four data points measure 
the divergence per unit time of the Tq — 8 calculation from 
(i) the quantum mechanics, (ii) the three-step semiclassical 
approximation, (iii) the five-step approximation, and (iv) the 
seven-step approximation. The upper and lower sets of points 
correspond to Heisenberg times N = 128 and -/V = 256, re- 
spectively. Theoretical curves corresponding to Eq. with 
c = 0.29, are drawn through the data. 

FIG. 3. Exponential convergence of the quantity 

Dt q = At q A^_ 1 towards the effective one-step long-time 
semiclassical propagator D* (Eq. ^). Data provided for 
Z>2 . . . -De ; Dg is used as the reference. 
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